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First-principles alloy theory, formulated within the exact muffin-tin orbitals method in combi¬ 
nation with the coherent-potential approximation, is used to study the mechanical properties of 
ferromagnetic body-centered cubic (bcc) Fei_i:M 2 , alloys (M=Mn or Ni, 0 < a: < 0.1). We consider 
several physical parameters accessible from ab initio calculations and their combinations in various 
phenomenological models to compare the effect of Mn and Ni on the properties of Fe. Alloying is 
found to slightly alter the lattice parameters and produce noticeable influence on elastic moduli. 

Both Mn and Ni decrease the surface energy and the unstable stacking fault energy associated with 
the {110} surface facet and the {110}(111) slip system, respectively. Nickel is found to produce 
larger effect on the planar fault energies than Mn. The semi-empirical ductility criteria by Rice 
and Pugh consistently predict that Ni enhances the ductility of Fe but give contradictory results in 
the case of Mn doping. The origin of the discrepancy between the two criteria is discussed and an 
alternative measure of the ductile-brittle behavior based on the theoretical cleavage strength and 
single-crystal shear modulus G{110}(111) is proposed. 

PACS numbers: 71.15.Nc, 62.20.-x, 71.20.Be, 75.50.Bb 


I. INTRODUCTION 

Duplex stainless steels (DSS) have achieved great suc¬ 
cess in the development of stainless steels and have led 
to numerous commercial applications. They consist of 
approximately equal amounts of ferrite and austenite. 
Though DSS successfully combined the merits of these 
two phases, the brittle cleavage problem in the ferrite 
component has attracted further consideration [l|. In 
order to grasp the mechanism of ductility in DSS, in the 
present work we concentrate on the effect of Mn and Ni 
additions on the intrinsic properties of ferromagnetic Fe. 
Both alloying elements are fundamental constituents of 
DSS and they are primarily responsible for the stabiliza¬ 
tion of the austenite phase in the low-carbon alloys. 

The mechanical properties of engineering alloys are de¬ 
termined by a series of mechanisms controlling the dislo¬ 
cation glide (such as elasticity, grain size, texture, grain 
boundary cohesion, solid solution and precipitation hard¬ 
ening). Describing such phenomena using merely ab ini¬ 
tio calculations is not feasible. On the other hand, the 
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mechanical properties are closely related to several mate¬ 
rial parameters which are within reach of such theoreti¬ 
cal tools. In particular, the fracture work depends on the 
surface energy, elastic energy, and plastic energy 0 . Fol¬ 
lowing the relationship between the crack and the surface 
energy (ts), and the plastic deformation and the unstable 
stacking fault (USF) energy ( 7 „), Rice Q proposed that 
the brittle-ductile behavior can be effectively classified 
by the Ts/Tu ratio. In body-centered cubic (bcc) system, 
plastic deformation shows various kinds of complexity 
due to its intricate dislocation core structure, multiple 
slip systems and the cross slip between them. According 
to the Peierls-Narbarro model, the shear stress necessary 
for moving a dislocation in the lattice can be estimated 
as Tp = 2Gj{l — ^) exp[—2a/(6(l — i^))], where G is the 
shear modulus, and v the Poisson ratio. Here a repre¬ 
sents the characteristic interlayer distance, and b is the 
Burgers vector. For a fixed structure, a/b depends only 
on the slip system. Accordingly, the greatest a and the 
shortest b give the lowest resistant stress, which corre¬ 
sponds to the most likely slip system. Usually, such slip 
systems have the closest packed atomic plane and the 
closest packed slip direction. In the bcc lattice, this is 
the {110}(111) slip system, i.e., slip occurs along the 
(111) direction on the {110} planes. A recent molecular 
dynamics study Q showed that other plane families such 
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as {112} and {123} can be also efficient in bcc Fe. 

Unfortunately, both the USF energy and the surface 
energy for a particular crystallographic facet are ex¬ 
tremely difficult to verify from experiments, and thus to¬ 
day the limited available data are exclusively based on 
first-principles quantum-mechanical calculations. How¬ 
ever, long before the era of “a 6 initio materials science”, 
Pugh summarized dozens of experimental data and de¬ 
veloped a simple phenomenological model 0 which cor¬ 
relates the fracture properties with some easily accessible 
elastic parameters. According to this model, the ductil¬ 
ity of materials can be characterized by the ratio between 
the bulk modulus B and the shear modulus G. Ductile 
materials usually possess B/G > yp = 1.75, whereas 
brittle materials have typically B/G < yp. The Pugh 
ratio is often used in modern computational materials 
science to assess and predict trends at relatively low cal¬ 
culation cost. 

One can easily establish a qualitative connection be¬ 
tween the Pugh criterion and that of Rice. The surface 
energy is known to correlate with the cohesive energy 
1 ^, which in turn is to a large degree proportional to the 
bulk modulus since larger cohesive energy requires larger 
curvature of the binding energy curve. The USF energy 
represents the energy barrier upon shifting one half of the 
crystal over the other half. The slope of the so defined 
7 -surface is expected to be proportional to (maximum 
of the curve) and also to the Peierls-Narbarro stress. As¬ 
suming a constant Poisson ratio, for a fixed lattice the 
latter is given merely by the shear modulus. Hence, at 
least on a qualitative level one may expect that "fs/"(u 
and B/G should follow similar trends, e.g., upon alloy¬ 
ing. The above correlations will be carefully investigated 
in the present work in the case of Fe-Mn and Fe-Ni alloys. 

For an isotropic polycrystalline material the two elas¬ 
tic moduli B and G determine the Poisson ratio v = 
{iB/G—2)/(&B/G+2). Consequently, the Pugh criterion 
imposes the constraint p > (3yp — 2)/(6yp -I- 2) ss 0.26 
for the Poisson ratio of ductile materials. The polycrys¬ 
talline bulk and shear moduli are derived from single¬ 
crystal elastic constants by Hill averaging of the Voigt 
and Reuss bounds 0 - In the case of cubic crystals, 
these two bounds are related to Pugh’s ratio by 

^ / 5 + P) ( 22-1 + 3)-i (Voigt) 

G \ i(|2-i+Ai)(3 + 22) (Reuss) ’ 

where 2 = 2Cii/{Gii — C 12 ) is the Zener anisotropy 
ratio. Hence, for both limiting cases, the Pugh ratio in¬ 
creases linearly with = Cx^/Caa,. For isotropic crystals 
2=1 and thus the Pugh criterion for ductility requires 
/i ^ 1.08. For weakly anisotropic lattices (2 ~ 2, which 
is the case of the present alloys), the above condition 
modifies to 77 > 0.92 or 77 > 1.07 corresponding to the 
Reuss or Voigt value, respectively. For this anisotropy 
level, the Hill average requires /r > 0.99 for ductile ma¬ 
terials. The above defined /r parameter is closely con¬ 
nected to the Cauchy pressure {C 12 — C 44 ) = C 44 (/i — 1). 
Negative Cauchy pressure (i.e., fj, < 1) implies angular 


character in the bonding and thus characterizes covalent 
crystals while positive Cauchy pressure (/r > 1) is typical 
for metallic crystals [HI- In other words, if the resistance 
of the crystal to shear (C 44 ) is larger than the resistance 
to volume change (<^ 12 ), the system shows covalent fea¬ 
tures and vice versa. This condition based on the Cauchy 
pressure is in line with the one formulated by Pugh. 

Compared to the criteria based on Pugh ratio and 
Cauchy pressure, the condition developed by Rice is ex¬ 
pected to give more physical insight into the mecha¬ 
nism of plastic deformation. However, the application 
of Rice condition based on planar defect energies is re¬ 
lated to the specific deformation mode and crystal ori¬ 
entation. In principle, one should analyze all possible 
cases with appropriate statistical weights. Except simple 
cases, such as the face-centered cubic (fee) lattice, this 
makes the Rice approach rather cumbersome. In turn, 
the Pugh rule was established from a series of experi¬ 
mental data and is based on simple physical quantities. 
Although many assumptions are involved in this empiri¬ 
cal approach, it has been confirmed to be practical in a 
large number of former scientific studies. 

In the present work, we focus on the alloying effect of 
Ni and Mn on ferrite Fe-based dilute alloys. In order 
to study the possible mechanisms behind the well-known 
facet cleavage problem in DSS alloys, which to a large ex¬ 
tent limits their engineering applications, here we carry 
out systematic electronic structure calculations to find 
the single-crystal and polycrystalline elastic parameters 
and the surface and USF energies as a function of chemi¬ 
cal composition. We consider the {110}(111) slip system 
due to its special role in the plastic deformation of the 
bcc lattice. Accordingly, we compute the surface energy 
for the {110} crystallographic facet and the USF energy 
for the {110} plane along the (111) slip direction. The 
theoretical predictions are used to assess various ductility 
parameters discussed above. 

The rest of the paper is divided into three main sec¬ 
tions and conclusions. Section [IT] gives a brief description 
of the employed theoretical tools, an introduction to the 
theory of elasticity and planar defects, and lists the im¬ 
portant numerical details. The results are presented and 
discussed in Sections Ell and El respectively. 


II. THEORETICAL TOOLS 
A. Total energy method 

The total energy calculations were performed u sing 
the exact muffin-tin orbitals (EMTO) method [l^ [T3 - 
El- This density functional theory (DET) solver is sim¬ 
ilar to the screened Korringa-Kohn-Rostoker method, 
where the exact Kohn-Sham potential is represented by 
large overlapping potential spheres. Inside these spheres 
the potential is spherically symmetric and constant be¬ 
tween the spheres. It was shown pm, 0 that us¬ 
ing overlapping spheres gives a better representation of 
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the full-potential as compared to the traditional muffin- 
tin or atomic-sphere approximations. Within the EMTO 
method, the compositional disorder is treated using the 
coherent-potential approximation (CPA) [13) [13 the 
total ener^ is computed via the full-charge density tech¬ 
nique [1,13 [I3- The EMTO-CPA method has been in¬ 
volved in many successful applications focusing on the 
thermo-physical properties of alloys and compounds [23- 
1 ^ . surface energies [Sll Ull and stacking-fault ener¬ 
gies [33 - l35j| . Here we would like to make a brief com¬ 
ment on CPA. Due to the single-site nature, CPA cannot 
account for the effect of local concentration fluctuations, 
short-range order effect and local lattice relaxation. This 
means that modeling the alloy by CPA one implicitly 
assumes a completely homogeneous random solid solu¬ 
tion with a rigid underlying crystal lattice. On the other 
hand, CPA can be used to describe fluctuations within 
the supercell, e.g. near the unstable stacking fault or 
surface, by using variable concentrations for each atomic 
layer. This feature is used here as well when studying the 
segregation effects (see Section II.C). In the case of Fe- 
Ni and Fe-Mn alloys, embedding a single impurity atom 
into the mean-field (CPA) Green function means that we 
omit the inter-dependence between magnetic state and 
the local environment, which in turn could lead to spe¬ 
cific local relaxation effects. In particular, within the 
present approximation the magnetic state of Mn is de¬ 
termined by the coherent Green function (representing 
the average effect of all alloying elements) rather than by 
the actual nearest neighbor atoms. Nevertheless, such 
local effects are expected to be relatively small when the 
impurity level is below ~ 10 % (i.e. all impurities are 
mostly surrounded by Fe). 

The self-consistent EMTO calculations were performed 
within the generalized-gradient approximation proposed 
by Perdew, Burke, and Ernzerhof (PBE) [s^]. The per¬ 
formance of this approximation was verified for the Fe- 
based systems in many former investigations [s^ - l^ . 
The one-electron equations were solved within the soft¬ 
core scheme and using the scalar-relativistic approxima¬ 
tion. The Green’s function was calculated for 16 com¬ 
plex energy points distributed exponentially on a semi¬ 
circular contour including states within 1 Ry below the 
Fermi level. In the basis set, we included s, p, d, and 
/ orbitals and for the one-center expansion of the full- 
charge density an orbital cut-off Imax = 8 was used. The 
total energy was evaluated by the shape function tech¬ 
nique with cut-off = 30 M- For the undistorted 

bcc structure, we found that a homogeneous fc-mesh of 
37 X 37 X 37 ensured the required accuracy. For the elas¬ 
tic constant calculations, we used about ~ 30000 — 31000 
uniformly distributed k points in the Brillouin zones of 
the orthorhombic and monoclinic structures. The po¬ 
tential sphere radii for the alloy components were fixed 
to the average Wigner-Seitz radius. For the surface en¬ 
ergy and USF energy calculations, the A:-mesh was set to 
9 X 23 X 2 in the the Brillouin zone of the base-centered 
orthorhombic super-cell. 


The impurity problem was solved within the single-site 
CPA approximation, and hence the Coulomb system of a 
particular alloy component i may contain a non-zero net 
charge. Here the effect of charge misfit was taken into ac¬ 
count using the screened impurity model (SIM) [43.1^. 
According to that, the additional shift in the one-electron 
potential and the corresponding correction to the total 
energy are controlled by the dimensionless screening pa¬ 
rameters ai and /3, respectively. The parameters ai are 
usually determined from the average net charges and elec¬ 
trostatic potentials of the alloy components obtained in 
regular supercell calculations [4l|. The second dimen¬ 
sionless parameter /3 is determined from the condition 
that the total energy calculated within the CPA should 
match the total energy of the alloy obtained using the su¬ 
percell technique. For most of the alloys, the suggested 
optimal values of ai and j5 are ^ 0.6 and ~ 1.2, respec¬ 
tively [13,1111. Often ai for different alloy components 
are assumed to the same and ,5 = 1 is used. 

In the present work, considering the volume sensitivity 
of ai, a dynamic SIM [131 was applied for the Fe-Ni 
system. According to that, the optimal ai varies as 
1.099,1.130,1.160,1.179,1.186,1.182,1.179,1.172,1.163 
for 9 lattice parameters a ranging from ^ 2.805 to 
^ 2.891 A (with interval of ~ 0.011 A). The second 
SIM parameter /3 was fixed to 1. For Fe-Mn, the SIM 
parameters show rather weak volume dependence, so 
ai and a = fdui were chosen to be 0.79 and 0.90, 
respectively. The above figures for the SIM parameters 
were obtained for a specific concentration (6.25 % Ni or 
Mn in Fe) and thus they are strictly valid only for that 
particular disordered system. However, in this study 
we assumed the same values for all binaries, i.e., the 
concentration dependence of the SIM parameters was 
omitted. 


B. Elastic properties 

The elastic properties of a single-crystal are described 
by the elements of the elasticity tensor. In a cubic lattice, 
there are three independent elastic constants: Cu, C 12 
and C 44 . The tetragonal shear elastic constant (C") and 
the bulk modulus {B) are connected to the single-crystal 
elastic constants through B = (Cu -I- 2 Ci 2)/3 and C' = 
(Cii - Ci2)/2. 

The adiabatic elastic constants are defined as the sec¬ 
ond order derivatives of the energy (E) with respect to 
strain. Accordingly, the most straightforward way to ob¬ 
tain the elastic parameters is to strain the lattice and 
evaluate the total energy change as a function of lattice 
distortion. In practice, the bulk modulus is derived from 
the equation of state, obtained by fitting the total energy 
data calculated for a series of volumes (V) or cubic lattice 
constants (a) by a Morse-type function [^. Since we are 
interested in equilibrium (zero pressure) elastic parame¬ 
ters, the fitting should involve lattice points nearby the 
equilibrium. It is known that bcc Fe undergoes a weak 
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FIG. 1: (Color online) Schematics of the model structures 
used to compute the USF energy (panel a) and the surface 
energy (panel b). See text for details. 


magnetic transition at slightly enlarged lattice parame¬ 
ters (a ~ 2.9 A), which should be excluded from the Ht- 
ting interval to avoid undesirable scatter of the computed 
bulk modulus Q. Accordingly, the present Morse-type 
ht was limited to lattice parameters smaller than 2.9 A. 
The same interval was applied to all considered binary 
alloys to minimize the numerical noise due to the volume 
mesh. 

Since the total energy depends on volume much more 
strongly than on small lattice strains, volume-conserving 
distortion are usually more appropriate to calculate C 
and Cii- Here, we employed the following orthorhombic 
{Do) and monoclinic {Dm) deformations 
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where <5 denotes the distortion parameter. For both 
lattice deformations, six different distortions 5 = 

0,0.01, 0.02,..., 0.05 were used. These deformations lead 
to total energy changes AE{6o) = 2VC'So + 0{5o) and 
AE{6m) = 2 VC44Sm + 0{6m), respectively, where O 
stands for the neglected terms. 

The polycrystalline shear modulus was obtained from 
the single-crystalline elastic constants according to the 
Hill average [i^, G = {Gy + Gii)/2, of the Voigt, 
Gv = (Cii — Gi 2 + 3 C' 44 )/ 5 , and Reuss, Gr = 5(C'ii — 
C\2)G44l\^G44 + 3(C'ii - C' 12 )], bounds. 


C. Surface and unstable stacking fault energy 

The schematics of the USF and surface model are illus¬ 
trated in Fig.[T] For USF (panel a), we show a top view 
along the [ 110 ] direction of the two atomic layers next 
to the stacking fault. The slip direction for the second 


(filled open circles) atomic layer is also shown. For the 
surface model (panel b), we show the side view along the 
[ 001 ] direction of the two surface atomic layers (filled cir¬ 
cles) and the two empty layers (empty circles) mimicking 
the vacuum. 

The surface energy ( 7 ^) of chemically homogeneous al¬ 
loys (without segregation) was extracted from total en¬ 
ergy calculations for slabs following the method proposed 
by Fiorentini and Methfessel [i^. No reference to sepa¬ 
rately computed bulk total energies is made in this ap¬ 
proach. Instead, the reference bulk energy per layer is 
deduced from a linear fit to the slab total energy ver¬ 
sus the slab thickness. Slabs with 8 , 10, and 12 atomic 
layers separated by 6 vacuum layers were used. Accord¬ 
ing to Punkkinen et al. [l^ the surface relaxation effects 
are small for the {110} surface facet of Fe and have an 
insignificant effect on the magnitude of the surface en¬ 
ergy. Our calculations confirmed this finding, and there¬ 
fore here we omitted the effect of relaxation in all surface 
energy calculations. 

The USF energy ( 7 ^) was obtained in the following 
way. First, the total energies of 16, 20, and 24 layers 
thick supercells (each containing two USF) were calcu¬ 
lated, and the bulk energy per layer was extracted from 
a linear fit to the energy versus thickness. Because the 
USF seriously disturbs the atomic arrangement in the 
vicinity of the stacking fault plane, structural relaxations 
perpendicular to the fault plane have to be taken into ac¬ 
count. We considered the relaxation effect for pure bcc 
Fe and found that the first inter-layer spacing at the USF 
increases by approximately 4%. For the Fe-based alloys, 
we kept the Mn/Ni content below 10%, which is expected 
to have a small additional effect on the layer relaxation. 
Therefore, we used the structural relaxation obtained for 
pure Fe for all binaries considered here. 

We also considered the effect of segregation on the sur¬ 
face energy and the USF energy. To this end, we varied 
the composition of the atomic layers next to the stack¬ 
ing fault or at the surface. Unfortunately, the above de¬ 
scribed method to compute the formation energies based 
entirely on the slabs/supercells with different sizes is 
not applicable when the chemical concentration profile 
is no longer homogeneous. Instead, all segregation stud¬ 
ies were performed for a slab/supercell with fixed size 
and the impact of segregation was interpreted with re¬ 
spect to the homogeneous slab/supercell with the same 
size. In that way, possible numerical errors associated 
with the reference system were kept at minimal level. In 
particular, the effect of segregation on the USF energy 
was derived by subtracting the energy of the supercell 
without USF from the energy of the supercell with USF, 
both with segregation on the slip planes. The surface en¬ 
ergy for the chemically inhomogeneous systems was com¬ 
puted following the method proposed in Ref. [i^ based 
on the grand canonical ensemble and involving the effec¬ 
tive chemical potentials. 
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III. RESULTS: INTRINSIC MATERIAL 
PARAMETERS 

A. Parameters of ferromagnetic bcc Fe 

We assess the accuracy of the EMTO method for the 
present systems by comparing the results obtained for 
pure Fe with those derived from previous full-potential 
calculations and the available experimental data. Re¬ 
sults are listed in Table ID The present lattice parame¬ 
ter a is consistent with those from Refs, [i^, [s^. All 
theoretically predicted equilibrium lattice parameter are, 
however, shghtly smaller than the experimental value, 
2.866 A [^, which is due to the weak PBE over binding. 
Theory reproduces the experimental magnetic moment 
accurately. 

In Table m the present single-crystal elastic constants 
and the bulk modulus of bcc Fe are compared with data 
from literature. The EMTO results are consistent with 
the other theoretical elastic constants, the deviations be¬ 
ing typical to the errors associated with such calculations. 
However, compared to the experimental data, EMTO is 
found to overestimate Cn by about ^ 19%, and under¬ 
estimates C 44 by ^ 18%. The larger Cn (and the larger 
bulk modulus) can be partly accounted for by the un¬ 
derestimated lattice parameter [s^. We notice that the 
excellent agreement in Ref. between the theoretical 
and experimental bulk moduli is somewhat suspicious 
since the quoted single-crystal elastic constants results 
in 187 GPa for the bulk modulus (compared to 172 GPa 
reported in Ref. [1^), which is close to the present fin ding 
and also to that from Ref. [1^ . The pseudopotential 
equilibrium lattice parameter is the largest among the 
theoretical predictions listed in the table, which might 
explain why the pseudopotential bulk modulus is sur¬ 
prisingly close to the experimental value. 

The EMTO surface energy is in close agreement with 
the full-potential result and also with the semi- 
empirical value valid for an “average” surface facet [s^. 
The present USE energy (1.08 J/m^) reproduces the full- 
potential result from Ref. within ~ 10%, which is 
well acceptable taking account that the present method 
is based on muffin-tin and full-charge density approxi¬ 
mations. We recall that today it is not yet possible to 
measure the USE energy. 


B. Equilibrium volumes and magnetic moments of 
Fe-Ni and Fe-Mn alloys 

According to the equilibrium phase diagrams, in fer¬ 
romagnetic bcc Fe the maximum solubility of Mn is 
^ 3% and that of Ni is 5.5%. Nevertheless, in du¬ 
plex stainless steels, solid solutions with substantially 
larger amount of Mn and Ni may appear as a result of 
the delicate balance between volume and surface effects. 
Hence, in the present theoretical study we go well be¬ 
yond the experimental solubility limits and discuss the 


trends up to 10% solute atoms in the Fe matrix. The 
calculated equilibrium lattice parameters for bcc binary 
Fei_a;Mna; and Fei_a;Nia; alloys are shown in Fig. [5] for 
0 < 2 : < 0.1. We find that for both binaries the lattice 
parameters show a weak dependence on the amount of 
alloying addition x. Our results, especially the alloying 
trends, are in reasonable agreement with the available 
experimental data [H, • 

The lattice parameter of bcc Fe increases with Mn ad¬ 
dition as X increases from 0 to 0.05. When Fe is alloyed 
with more than 5 % Mn, the lattice parameter shows a 
small negative slope as a function of x. Due to the lim¬ 
ited solubility of Mn in bcc Fe, the experimental data is 
restricted below ~ 3% Mn. In contrast to Mn, the lat¬ 
tice parameter of Fei_a;Nia; monotonically increases with 
the Ni amount. Sutton et al. reported that the lat¬ 
tice parameter of Fei-^Nij, is a linear function of x up 
to 5% Ni with small slope and remains constant at 
larger concentrations. This trend is reproduced by the 
present theory, although the switch from a weak positive 
slope to a constant value seems to be shifted to larger 
concentrations (Fig. [2]). 

The present theory slightly overestimates the observed 
volume expansion due to alloying. This can have both 
theoretical and experimental origin. From the side of 
theory, the exchange-correlation functional (PBE) may 
produce a noticeable effect since it performs differently 
for the present metals. For instance, while PBE under¬ 
estimates the equilibrium lattice parameter of bcc Fe by 
^ 1.0%, it slightly overestimates that of fee Ni ('^ 0.3%) 
[13, [ 5 ^. This error results in ^ 0.004 A “extra” increase 
of the lattice parameter at 10% Ni compared to the ideal 
case when the DFT error stays at the same level for dif¬ 
ferent elements. A careful inspection of Fig. [^shows that 
the above error accounts for about 30% of the deviation 



FIG. 2: (Color online) Lattice parameters (in units of A) of 
bcc Fei_i:Mn 2 : and Fei-^NL alloys as a function of compo¬ 
sition. For comparison, the measured values by Sutton 
and Owen [s^] are also included. 
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TABLE I: Theoretical and experimental lattice parameter a (in units of A), magnetic moment /r (in /is), single-crystal elastic 
constants (in GPa), surface energy 7 s (in J/m^) and USE energy 7 ^ (in J/m^) for ferromagnetic bcc Fe. The quoted theoretical 
methods are EMTO: present results; PAW: full-potential projector augmented wave method; PP: ultrasoft non-norm-conserving 
pseudopotential; FLAPW: all-electron full-potential linear augmented plane wave method; FPLMTO: full-potential linear 
muffin-tin orbitals method. The low-temperature (4 K) experimental elastic constants are listed for comparison, as well as the 
semi-empirical surface energy. 


Method 

a 


Cii 

C 12 

C44 

C 

B 

Is 

7 u 

EMTO 

2.838 

2.21 

288.9 

133.7 

100.0 

77.6 

185.4 

2.47 

1.08 

PAW 

2.838[4^ 

2 . 21 [^ 

271 [4^ 

145[^ 

101 [^ 

63 

172 [4^ 

2.50[T7] 

0.98[^ 

PP 

2.848[^ 

2.25[^ 

289 

118[^ 

115[^ 

85.5[^ 

176 



FLAPW 

2.84[^ 

2.17[^ 

279 

140[^ 

99 

69 

186[^ 



FPLMTO 

2.812[M1 


303 

150[^ 

126[^ 

76.5[^ 

201 [^ 



Expt. 

2 .866 [511 

2.22[511 

243.1[^ 

138.1[511 

121.9[511 

52.5[^ 

173.1[^ 

2.41 



obtained between the theoretical and experimental slopes 
of a{x) for the Fe-Ni system. Another possible source of 
error is the neglect of all thermal effects in the present 
theory. From the experimental point of view, not all mea¬ 
sured data points are relevant for the present comparison. 
That is because both Ni and Mn are fee stabilizer, and 
thus a perfect bcc single-phase is hard to attain. More¬ 
over, in the process of preparing the samples, different 
quenching temperatures and various annealing methods 
lead to subtle changes in the measured lattice parameter. 

The overall small influence of Mn/Ni alloying addition 
on the lattice parameter of Fe can be explained from the 
atomic volume and the magnetic moment. The atomic 
volumes of Mn and Ni are very close to that of Fe. Hence, 
based on the linear rule of mixing, small amounts of Mn 
and Ni should not affect the volume of Fe to a large ex¬ 
tent. On the other hand, the observed small increment 
can be ascribed to the complex magnetic interaction be¬ 
tween the solute atoms and ferromagnetic Fe matrix. 

The local magnetic moments in Fe-Mn and Fe-Ni are 
displayed in Fig. |3]as a function of composition and lat¬ 
tice parameter. These results were extracted from the 
CPA calculations and the local magnetic moments rep¬ 
resent the magnetic moment density integrated within 
the corresponding Wigner-Seitz cells. Around the equi¬ 
librium lattice parameters (~ 2.84 — 2.85 A), the Mn 
moments are antiparallel to those of Fe. The magni¬ 
tude of Mn local magnetic moment ^(Mn) decreases from 
^ 1.7/rB to about zero as the amount of Mn increases 
from zero to 10% (Fig. [31 upper right panel). This 
trend is consistent with the previous findings, namely 
that above ^ 10% Mn, the coupling between Fe and Mn 
becomes ferromagnetic [s^. At the same time, the mag¬ 
netic moment of Fe /i(Fe) is increased slightly with Mn 
addition up to ^ 5% Mn, above which a very weak de¬ 
crease of At(Fe) can be observed (Fig. [31 upper left panel). 
The enhanced /r(Fe) expands the volume as a result of 
the positive excess magnetic pressure [o^- For x > 0.05, 
the small negative slope of /i(Fe) leads to vanishing ex¬ 
cess magnetic pressure and thus to shrinking volume, in 
line with Fig. [3] On the other hand, the Ni moments 
/r(Ni) always couple parallel with the Fe moments (Fig. 



FIG. 3: (Color online) The map of the local magnetic mo¬ 
ments fi{x) (in units of hb) as a function of chemical com¬ 
position and lattice parameter a (in units of A) for Fe and 
Mn in Fei_a;Mna; (upper panels) and Fe and Ni in Fei-^Nb 
(lower panels). 


m lower right panel) and Ni addition increases /r(Fe) (Fig. 
[31 lower left panel). This in turn increases the excess 
magnetic pressure and yields monotonously expanding 
volume with x. The electronic structure origin of the 
alloying-induced enhancement for the Fe magnetic mo¬ 
ment for the present binaries is discussed in Ref. [blj. 

In order to gain more insight into the magnetic cou¬ 
pling between Mn and the bcc Fe host, we constructed a 
2 x2x2 supercell in terms of the conventional bcc unit cell 
(16-atom supercell) containing one or two Mn impurity 
atoms, which correspond to 6.25% Mn and 12.5% Mn, re¬ 
spectively. In the case of one isolated impurity atom, Mn 
couples antiferromagnetically to Fe and the average mag¬ 
netic moments of Fe and Mn are almost identical to those 
from the corresponding CPA calculation. In the case of 
12.5% Mn in the supercell, both Mn moments align par¬ 
allel with respect to the surrounding Fe host irrespective 
of their relative positions in the supercell. However, the 
magnitude of their moments scales with the distance be- 
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tween them. If the two Mn atoms are located at nearest 
neighbor sites, their local magnetic moments are strongly 
reduced which may be interpreted as a result of compet¬ 
ing magnetic interactions between Mn-Mn and Mn-Fe. 
At the same time, the local magnetic moments of the 
nearest Fe atoms remain almost at the level of pure bcc 
Fe, similar to the CPA result. These tests, based on su¬ 
percell calculations, confirm that CPA accurately mimics 
the behavior of the magnetic moments versus concentra¬ 
tion in the dilute limit. 


C. Elastic parameters 

The theoretical single-crystal elastic constants and 
polycrystalline elastic moduli of the Fe-Ni and Fe-Mn 
alloys are plotted in Figs. 2] and [S] as a function of com¬ 
position and the corresponding data are listed in Tablellll 
The results indicate rather impressive alloying effects of 
Mn/Ni on the elastic parameters. We find that Mn and 
Ni produce similar effects on Cn of Fe. Namely, Cn de¬ 
creases with X below ^ 7% alloying addition, and then 
increases at larger concentrations. But Mn and Ni give 
different effects on Ci 2 . Below ^ 6% Ni, C 12 keeps con¬ 
stant with Ni content, and then strongly increases with 
further Ni addition. On the other hand, C 12 drops from 
133.7 to 103.2 GPa as the Mn concentration increases 
from zero to ~ 8%. At larger concentrations, the effect of 
Mn on C 12 changes sign. The peculiar concentration de¬ 
pendencies obeyed by Cn and C 12 originate mainly from 
the non-linear trend of the corresponding bulk modulus 
shown in Fig. O 

The two cubic shear elastic constants also exhibit com¬ 
plex trends. C decreases with Ni addition, which means 
that Ni decreases the mechanical stability of the bcc lat¬ 
tice. At the same time, Mn is found to have very small 
impact on the tetragonal shear elastic constant. Cn in¬ 
creases with Mn addition, but remains nearly constant 
upon alloying Fe with Ni. The present trends for the elas¬ 
tic constants are very close to those reported by Zhang 
et al. using the same total energy method. The 

small differences seen at larger concentrations, especially 
in the case of the bulk modulus, are due to the different 
numerical parameters used in these two works. 

The theoretical polycrystalline elastic moduli B, E, 
and G (Fig.lbl reproduce reasonably well the experimen¬ 
tal trends (62l [ except perhaps for alloys encompassing 
about 10% solute. This discrepancy is most likely due 
to the fact that the Fe-Mn and Fe-Ni binary alloys with 
X ^ 0.05 exist as a mixture of bcc and fee phases. On the 
theory side, the lattice parameters of Fe-Ni and Fe-Mn 
alloys are underestimated as a result of the employed 
exchange-correlation approximation (Fig. [5]), which at 
least partly explains why theory in general overestimates 
the elastic moduli. 

On a qualitative level, the trends of the bulk mod¬ 
uli we understand as the result of the interplay between 
chemical and volume effects. The calculated bulk mod- 
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FIG. 4: (Color online) Single-crystal elastic constants of bcc 
Fei-ajMna; and Fei-^Nh alloys as a function of composition. 



FIG. 5: (Color online) Poly crystalline elastic moduli of bcc 
Fei-ajMna; and Fei-^Nh alloys as a function of composition. 
The experimental data are from Ref. [^ . 


uli of pure Ni and Mn in the bcc lattice (i3(Mn)Ri 222 
GPa, B(Ni)Ki 193 GPa) are both larger than that of Fe. 
Hence, at large concentrations the bulk moduli of binary 
alloys should eventually increase. On the other hand, 
at low concentrations (less than ~ 7%), the interaction 
between the impurity atoms is weak and thus the trend 
of the bulk modulus is mainly governed by the volume 
effect. Increasing volume in turn produces a drop in the 
bulk modulus, in agreement with Fig. [5] 

The results for Fe-Ni indicate that the polycrystalline 
elastic moduli E and G monotonically decrease with so¬ 
lute concentration. For Fe-Mn, the alloying effect can be 
divided into two parts: when x is less than ~ 0.05, G 
and E remain constant with x. When x is larger than 
^ 0.05, both G and E increase with Mn content. Since 
the Voigt and Reuss bounds depend only on the single- 
















crystal shear elastic constants, the trends of G in Fig. 
[S] directly emerge from those of (744 and C (Fig. H]). 
We recall that for isotropic crystals, G reduces exactly 
to the single-crystal shear elastic constant {G^^ = G'). 
The Young modulus is a mixture of B and the Pugh 
ratio B/G, viz. E = ^B/{iB/G + 1), which is nicely 
reflected by the trends in Fig. [S] For Fei_a,Nia; with 
Ni content below ^ 6 %, B/G remains nearly constant 
with X (Table nil), and thus the corresponding E resem¬ 
bles the bulk modulus. At larger Ni concentrations, B/G 
increases substantially which explain the continuous de¬ 
crease of E. In the case of Fe-Mn, B/G decreases with 
Mn addition which gives a strong positive slope to the 
Young modulus as compared to that of B. 


D. Surface energy and unstable stacking fanlt 
energy 

The calculated surface energies and USF energies of 
Fei_a,Mn 2 , and Fei-^Nij, are shown in Fig. | 6 ]as a func¬ 
tion of X. Both Mn and Ni are predicted to decrease 
the surface energy of the {110} facet of bcc Fe. Nickel 
has a stronger effect than Mn. Namely, 10% Ni addition 
reduces the surface energy of Fe by 0.17 J/m^, which is 
about 7% of the surface energy of pure Fe, whereas Mn 
lowers the surface energy by 0.07 J/m^ (about 3%). The 
alloying effect on the surface energy can be understood 
on a qualitative level using the surface energies calculated 
for the alloy components [d^. Nickel has substantially 
smaller surface energy (considering the close-packed fee 
surfaces) than bcc Fe, and thus Ni addition is expected 
to reduce 7 s of Fe. The surface energy of a—Mn is larger 
than that of Fe However, when considering the bcc 
lattice, the surface energy of Mn turns out to be interme- 



X at. % Ni or Mn 


FIG. 6: (Color online) Surface energy and unstable stacking 
fault energy of bcc Fei_a;Mna; and Fei_a;Nb alloys as a func¬ 
tion of composition. The surface energy corresponds to the 
{110} facet and the USF energy to the {110}(111) slip system. 



X {at. % Mn) X {at, % Ni) 

FIG. 7: (Color online) Relative change ( —[7(a:^,3/) — 
7(0,0)]/7(0, 0)) of the USF energy and the surface energy 
as a function of bulk concentration (a;) of Mn/Ni (abscissae) 
and the additional solute concentration (y) at the fault plane 
(ordinates) for Fei-iMn^, and Fei_a;Ni3; (0 < x < 5%) alloys. 


diate between those of Ni and Fe, which is nicely reflected 
by the relative effects of Mn and Ni on 7 s. 

The USF energy shows a similar dependence on Ni/Mn 
alloying as the surface energy; it decreases from 1.08 to 
0.92 J/m^ upon 10% Ni addition, which is about 15% 
of the USF energy of pure Fe. Compared to the effect 
of Ni, 10% Mn produces a smaller change in 7 ^ (5.6%). 
Using the same methodology as for Fe and Fe-alloy, we 
computed the USF of hypothetical bcc Mn and bcc Ni for 
the present slip system and at the volume of bcc Fe. We 
obtained 0.63 J/m^ for Mn and 0.44 J/m^ for Ni. These 
figures are in line with the trends from Fig. [5] 

The segregation effect on the surface energy and USF 
energy are shown in Fig. [71 We find that both Mn and Ni 
segregate to the fault plane. Above ^ 1% solute concen¬ 
tration in bulk Fe, the surface energy and the USF energy 
decrease slightly with segregation of Ni/Mn to the layer 
at the planar fault. In this segregation calculation we 
kept the volume the same as in the bulk (host alloy). 
Therefore, only the chemical segregation effect is consid¬ 
ered and the local volume expansion effect nearby fault 
plane due to the segregating solute atoms was ignored. 
The results show that the segregation effect on the sur¬ 
face energy is larger for Ni than for Mn. Furthermore, 
the surface segregation effect slightly increases with in¬ 
creasing Mn content in Fe-Mn dilute alloy. For the USF 
energy, the segregation effect of Mn is similar to that 
of Ni when the concentrations are small. However, with 
increasing solution amount x, the segregation of Mn to 
stacking fault increases with high degree, but the effect 
of Ni remains almost the same. 

The surface energy and the USF energy are primarily 
determined by the properties of the surface layer and the 
slip layers, respectively. Segregation changes the concen¬ 
tration of the solute at the surface or slip plane. Because 
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TABLE II: Theoretical single-crystal elastic constants < 711 ( 0 :), Ci 2 {x), C 4 ,a{x), C'{x), polycrystalline bulk B(x), shear G{x) 
and Young E{x) moduli (GPa), Poisson ratio v, Pugh ductile/brittle ratio BjG, Cauchy pressure {C12 — G44) (GPa) and Zener 
anisotropy ratio G44/G' for Fei_a:Mna; (upper panel) and Fei-ajNia, (lower panel) alloys as a function of composition. 


X 

B 

G' 

Gii 

C12 

<^•44 

G 

E 

V 

(C42 - G44) 

B/G 

G44/C' 

FB 1 a; A/[ n j; 

0 

185.4 

77.6 

288.9 

133.7 

100.0 

90.3 

233.1 

0.290 

33.71 

2.05 

1.29 

0.01 

179.1 

76.0 

280.5 

128.4 

100.5 

89.9 

231.0 

0.285 

27.94 

1.99 

1.32 

0.02 

172.0 

75.2 

272.3 

121.8 

102.8 

90.7 

231.4 

0.276 

19.06 

1.90 

1.37 

0.03 

166.8 

74.1 

265.6 

117.4 

105.4 

91.5 

232.1 

0.268 

12.01 

1.82 

1.42 

0.04 

161.0 

73.3 

258.8 

112.1 

107.7 

92.3 

232.5 

0.259 

4.41 

1.74 

1.47 

0.05 

159.7 

73.1 

257.2 

111.0 

110.4 

93.6 

234.8 

0.255 

0.59 

1.71 

1.51 

0.06 

155.4 

73.6 

253.5 

106.3 

113.0 

95.2 

237.1 

0.246 

-6.71 

1.63 

1.54 

0.07 

153.1 

74.6 

252.5 

103.4 

115.1 

96.7 

239.6 

0.239 

-11.68 

1.58 

1.54 

0.08 

153.6 

75.6 

254.4 

103.2 

116.8 

98.1 

242.6 

0.237 

-13.55 

1.57 

1.54 

0.09 

157.7 

76.3 

259.4 

106.8 

117.9 

99.0 

245.7 

0.240 

-11.07 

1.59 

1.55 

0.1 

163.5 

76.8 

265.9 

112.3 

119.0 

99.8 

248.9 

0.246 

-6.71 

1.64 

1.55 

F6i — xNia; 

0 

185.4 

77.6 

288.9 

133.7 

100.0 

90.3 

233.1 

0.290 

33.7 

2.05 

1.29 

0.01 

183.6 

74.6 

283.0 

133.9 

97.3 

87.5 

226.5 

0.294 

36.6 

2.10 

1.31 

0.02 

181.2 

72.6 

278.0 

132.8 

96.6 

86.1 

223.1 

0.295 

36.2 

2.10 

1.33 

0.03 

177.7 

70.6 

271.8 

130.6 

96.3 

85.0 

220.0 

0.294 

34.4 

2.09 

1.36 

0.04 

174.6 

68.7 

266.1 

128.8 

96.4 

84.2 

217.5 

0.292 

32.4 

2.07 

1.40 

0.05 

173.3 

66.8 

262.4 

128.8 

97.0 

83.5 

215.9 

0.292 

31.7 

2.07 

1.45 

0.06 

173.6 

64.9 

260.2 

130.3 

98.2 

83.2 

215.2 

0.293 

32.1 

2.09 

1.51 

0.07 

176.8 

62.9 

260.6 

134.9 

99.4 

82.7 

214.6 

0.298 

35.5 

2.14 

1.58 

0.08 

181.8 

60.7 

262.8 

141.3 

100.6 

82.2 

214.2 

0.304 

40.7 

2.21 

1.66 

0.09 

188.3 

58.6 

266.4 

149.3 

102.2 

81.7 

214.2 

0.310 

47.1 

2.30 

1.74 

0.1 

196.9 

56.4 

272.1 

159.3 

103.4 

81.1 

213.9 

0.319 

55.9 

2.43 

1.83 


Mn and Ni have lower planar fault energies than Fe (in 
the bcc structure), both formation energies of the chem¬ 
ically homogeneous Fe-alloys are expected to decrease as 
the Mn or Ni solute concentration increases, in line with 
the present results. The effect is more pronounced for Ni 
since its USF/surface energy is smaller than that of Mn. 


IV. DISCUSSION 
A. Ductile and brittle properties 

The fracture behavior of a specific material rests on 
various conditions, such as, the presence of flaws, the 
way and the magnitude of the applied stress, tempera¬ 
ture, strain rate, and alloying elements. Alloying changes 
the interaction among atoms, and produces indispensable 
influences on the mechanical properties like the elastic 
response or the ductile versus brittle behavior. In the 
following, the ductility of the dilute Fe-Mn and Fe-Ni 
binary alloys is addressed using previously established 
effective theoretical models and phenomenological rela¬ 
tionships based on planar fault energies and elastic con¬ 
stants. These models are widely used and often overused 
in combination with ab initio calculations and thus a 
careful assessment of their performance is of common 
interest. Here we employ four criteria discussed in the 
introduction to estimate the effect of Mn/Ni on the duc¬ 


tility of ferromagnetic bcc Fe. These criteria are based on 
the Poisson ratio {v), the Cauchy pressure {C12 — C44), 
the Pugh ratio B/G, and the Rice ratio 7s/T u- We re¬ 
call that the first two criteria are closely connected with 
the Pugh conditions and thus no substantial deviations 
between them are expected. 

The present theoretical predictions are shown in Fig.|51 
Both B/G and Poisson’s ratio indicate that Mn makes 
the ferrite Fe-based alloys more brittle. According to 
Pugh criterion, about 4% Mn is needed to transfer the 
Fe alloy from the ductile into the brittle regime. On the 
other hand, small Ni addition (< 6%) keeps the good 
ductility of Fe whereas larger amounts of Ni make the 
Fe-Ni system more ductile. 

The Cauchy pressure {G12 — G44) follows very closely 
the trend of B/G and that of the Poisson ratio. How¬ 
ever, the Cauchy pressure becomes negative at a slightly 
larger concentration (~ 6%) than the critical value in 
terms of B/G (^ 4%). The small difference between the 
’’predictions” based on these phenomenological correla¬ 
tions originates from the elastic anisotropy of the present 
alloys (cf. Eq. (P)). 

According to the theory developed by Rice Q , the ratio 
Is/lu is associated with the fracture behavior. The in¬ 
crement of ^s/lu stimulates the creation of dislocations. 
In such cases, the stress around a crack tip will be re¬ 
leased by slipping the atomic layers. A decreasing ratio, 
on the other hand, means that the material will crack by 
opening new micro surfaces. 
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FIG. 8: (Color online) Poisson ratio v, Cauchy pressure 
(C \2 — C44) (in units of GPa), Pugh ratio B/G, and Rice 
ratio ')sl^u of bcc Fei-xMnj; and Fei_a;Nia; alloys as a func¬ 
tion of composition. 


Figure |8] shows that for Fei_a:Nia;, "fs/lu increases as 
a function of x indicating that the system becomes more 
ductile with Ni addition. Although Ni decreases both 
7 s and 7 k simultaneously, the alloying effect on the USF 
energy is more pronounced, resulting in an increasing 
Rice ratio from 2.27 corresponding to pure Fe to 2.47 
belonging to Feo.gNio.i- In the case of Mn addition, we 
observe relatively small changes in the Rice ratio. Small 
amounts of Mn make the Fe-Mn alloy slightly more duc¬ 
tile in terms of the Rice parameter (compared to pure 
Fe), but with increasing Mn content beyond ^ 6% the 
Ishu ratio saturates around 2.35. 

According to the present theoretical results for the Rice 
parameter, Ni has stronger effect in making the alloy duc¬ 
tile compared to Mn. In addition, more than 6% addition 
of Mn makes the bcc phase relatively more brittle. In 
fact, in terms of the USF, one would predict that both 
Ni and Mn make the dislocation formation in Fe more 
easy, and the effect of Ni is superior to that of Mn. On 
the other hand, both element decrease the surface energy 
as well, making the crack opening more likely. 

Before turning to the comparison between the Pugh 
and Rice criteria, we make an observation based on the 
present segregation studies. Allowing for surface and in¬ 
terface segregation leads to small changes in the USF en¬ 
ergy and the surface energy (cf. Section III.D). In terms 
of the Rice parameter, in alloys containing 5% impurity 
5% surface/USF segregation increases slightly the "fs/lu 
ratio (^ 2% for Fe-Ni and ~ 4% for Fe-Mn). That is 
because in both binary systems the USF energy is low¬ 
ered by a larger degree than the surface energy upon 
segregation (Fig. [T]). These changes are far below those 
associated with the effect of bulk concentration (Fig. [8]). 
In addition, taking into account the different time-scales 
for atomic diffusion and dislocation movement, we con¬ 
clude that the segregation effects may safely be omitted 


for the present discussion. 


B. Pugh criterion versus Rice criterion 

In the case of Fe-Ni alloys, the conclusion drawn from 
la/lu is consistent with the other three criteria based on 
R/G, Poisson ratio, and Cauchy pressure. When we look 
for the individual trends, we find that this consistency is 
to some extent incidental, especially at large Ni content. 
Namely, while both G and decrease with Ni addition, 
the trends for the surface energy and B strongly devi¬ 
ate from each other above ~ 5% Ni. At low Ni levels B 
and 7 s follow similar trends, but the relatively large B of 
Feo.gNio.i is not supported by the strongly decreased 7 s 
calculated for this alloy. The reason why 7s/7« and B/G 
still predict similar effects (from the Rice and Pugh con¬ 
ditions) is simply due to the strong decrease calculated 
for 7 k upon Ni addition. 

In the case of Fei-aiMna,, the conclusions made based 
on the Rice and Pugh criteria contradict each other. 
While B/G decreases for x < 0.08, "fs/lu shows a weak 
increase for these alloys. We find that for this system, G 
and 7 k follow completely different trends as a function of 
Mn content (Figs. [Slandin]), and the deviation between 
the trends of 7 s and B is also pronounced (although at 
low and intermediate Mn level both of them show de¬ 
creasing trends). From these results, we conclude that 
the two ductility criteria (Pugh and Rice) lead to incon¬ 
sistent results, a fact that questions their reliability and 
limits their scope. 

In the following we make an attempt to understand the 
origin of this discrepancy. To this end, we make use of 
the particular shear elastic constant associated with the 
present slip system as well as of the concept of theoretical 
cleavage stress. These two quantities are considered here 
as possible alternative measures of the materials resis¬ 
tance to dislocation slip and cleavage, respectively, and 
are expected to give a better estimate of the correspond¬ 
ing effects compared to the polycrystalline G and B em¬ 
ployed in the Pugh criterion. 

Within the Griffith theory of brittle fracture, the the¬ 
oretical cleavage stress is often approximated as 


rj -i f Biyjiji'y I mn N ,c\\ 

aci.{lmn\ = 1 — - I , ( 2 ) 

\ ^Imn J 

where {Imn} stands for the cleavage plane, Eimn, limn 
and dimn are the corresponding Young modulus, sur¬ 
face energy and interlayer distance, respectively. Ir- 
wing and Orowan extended the above equation (mod¬ 
ified Griffith equation) by including the plastic work 
before fracture. Here we neglect this additional term, 
*•6. limn represents merely the solid vacuum inter¬ 
face energy. Using our calculated elastic constants, lat¬ 
tice parameters and surface energy, we computed ad. 
for Fe-Mn and Fe-Ni for the {110} plane, for which 
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FIG. 9: (Color online) Relative changes {r; in %) for the sur¬ 
face energy and the theoretical cleavage stress for the {110} 
surface (upper panel) and for the USF energy and the single¬ 
crystal shear elastic modulus associated with the 1/2(111) slip 
system (lower panel) for bcc Fei-irMn^: and Fei-^Nia; alloys 
as a function of composition. Fixed-spin results are shown for 
the surface energy and the USF energy of the Fe-Mn system. 


Eno = 12C'C44B/(C'iiC44 + 3C'B). The alloying- 
induced changes for crci.jllO} are compared to those cal¬ 
culated for the surface energy in Fig.O upper panel. Here 
the change rjx{x) = [X{x) — X(0)]/X(0) {X{x) stands 
for the cleavage stress or the surface energy for Fei_a;Ma; 
alloy) is expressed relative to the corresponding value in 
pure Fe X(0). It is found that (x) and follow 
similar trends for Fe-Ni, but strongly deviate for Fe-Mn. 
Before explaining this deviation in the case of Mn dop¬ 
ing, we introduce the shear modulus associated with the 
present slip system. 

In bcc alloys, slip occurs primarily in the {110} plane 
along the (111) with Burgers vector (1/2,1/2,1/2). The 
associated shear modulus can be expressed as 


Fe-Mn Fe-Ni 



FIG. 10: (Color online) Local magnetic moments of Fe (upper 
curves) and Mn/Ni (lower curves) for Fei-^Muj, (left panels) 
and Fei-a;Ni 3 ; (right panels) as a function of layer index. Re¬ 
sults are shown for the model systems used for the USF energy 
(upper panels) and surface energy (lower panels) calculations. 
Different symbols correspond to various impurity levels (a:) as 
shown in the legend. 


We note that the above modulus in fact expresses 
the shear for any possible shear plane {Imn} which 
contains the (111) shear direction. In the original 
Pugh criteria, the averaged shear modulus G is used 
which in anisotropic materials can substantially differ 
from G'{110}(111). The alloying-induced changes for 
G'{110}(111) are compared to those calculated for the 
USF energy in Fig. [HI lower panel. Interestingly, we 
find that ryG{iio}<iii> and are practically identical 
for Fe-Ni. Both G{110}(111) and 7 „ decrease by ^ 15% 
when 10% Ni is added to Fe. Nickel substantially softens 
the elastic modulus associated with the {110}(111) shear 
which is nicely reflected by the decrease of the energy 
barrier for the {110}(111) slip. For Fe-Mn, ?7g{iio}(iii) 
and are also close to each other for Mn content below 
^5%, but show large deviations at larger concentrations. 
Adding more Mn to Feo.gsMno.os further decreases the 
USF energy, but G{110}(111) changes slope and shows 
a weak increase for Feo.gMno.i relative to the value for 
pure Fe. 

We suggest that the different behaviors obtained for 
Fe-Ni and Fe-Mn has to a large extent magnetic origin. 
That is because in contrast to Ni, Mn is a weakly itiner¬ 
ant magnet and thus any (here structural) perturbation 
can have a marked impact on its magnetic state. To il¬ 
lustrate that, in Fig. [TU] we plotted the local magnetic 
moments for the unit cells used for the USF energy and 
surface energy calculations. In the case of Fe-Ni, we see 
no substantial deviation in the local magnetic moments 
as we approach the planar fault area. Both Ni and Fe 
moments near the fault plane remain close to the bulk 
value. However, for Fe-Mn, the local magnetic moments 
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of Mn next to the planar fault prefer a very strong an¬ 
tiferromagnetic coupling with the Fe matrix irrespective 
of the bulk concentration. Although the bulk moments 
approach zero as the Mn content increases to 10%, the 
interface Mn moments remain around — 2.2/iB and those 
on the surface around —3.1/re- This strong antiferromag¬ 
netic Fe-Mn coupling near the USF interface and surface 
indicates an energetically stable configuration that can 
lower the corresponding formation energy. Indeed, re¬ 
moving this degree of freedom by fixing all Mn moments 
to the corresponding bulk value {i.e., modeling a situ¬ 
ation similar to the case of Fe-Ni system) increases the 
the surface energy and USF energy of Fe-Mn. The rel¬ 
ative changes of the corresponding yf® and 7 ^® values 
(FS stands for fixed-spin) relative to those of pure Fe are 
shown in Fig. [S] It is found that 7 ^® for Feo.gMno.i ap¬ 
proaches the USF energy of pure Fe. In fact, constraining 
the magnetic moments near the planar fault brings the 
trend followed by 7 ^® rather close to that of G{110}(111) 
(Fig. ini lower panel). A similar impact of the constrained 
magnetic moment is seen for the surface energy as well. 
Namely, for Fe-Mn yf ® and the theoretical cleavage stress 
show similar concentration dependencies (Fig. [HI upper 
panel). On this ground, we conclude that the deviations 
seen between the surface energy and cleavage stress and 
between the USF energy and the shear modulus in the 
case of Fe-Mn are due to the stable antiferromagnetic 
state of Mn near the planar faults. 

Using the FS results for the Rice parameter, we get 
a weak decrease of yf^/y^® above ^ 5% Mn addition 
to Fe (shown in Fig. [TTl upper panel). Thus, the fixed- 
moment result for the Rice parameter is somewhat closer 
to the original Pugh ratio, both of them predicting en¬ 
hanced intrinsic brittleness for the Fe-Mn solid solution. 
This demonstrates that the discrepancy seen between the 
two criteria (Fig. IH]) is partly due to the weak magnetic 
behavior of Mn associated with the planar defects but 
absent in the elastic moduli. We suspect that the situa¬ 
tion in non-magnetic alloys could be closer to the case of 
Fe-Ni, where a good parallelism between the Pugh and 
Rice conditions is found. However, further theoretical re¬ 
search is needed before making the final verdict for this 
question. 

Finally, we test the ratio between the cleavage energy 
and the shear modulus associated with the slip system, 
X{lmn} = aci.{lmn}/G{llQ}{lll), as a possible alterna¬ 
tive measure of the ductile-brittle behavior. In the upper 
panel of Fig. 1111 we compare the Rice parameters (from 
Fig. [5]) to A{110}. For Fe-Ni, we find a good correspon¬ 
dence between the two measures. Namely both of them 
increase with Ni addition, indicating enhanced ductility. 
For Fe-Mn, slightly larger deviation occurs at large Mn 
content. As we discussed above, part of this deviation 
may be ascribed to magnetism and in particular to the 
stable antiferromagnetic state of Mn around the planar 
faults. 

Since in bcc metals, cleavage predominantly occurs in 
the { 100 } plane, in addition to the previously discussed 
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FIG. 11: (Color online) Theoretical Rice parameter and the 
ratio between the theoretical cleavage stress and the shear 
modulus associated with the slip system (A) for bcc Fei_a:Mn 3 , 
and Fei_a:Nia; alloys as a function of composition. Fixed-spin 
results are shown for Rice parameter of the Fe-Mn system. 
Upper panel represents {110} as the cleavage plane, and lower 
panel {100} as the cleavage plane. 

cleavage stress, we also computed crcj,{ 100 } using the cor¬ 
responding surface energy and Uioo = dtCB/{Cn+Ci^)■ 
The so obtained A{100} parameters are shown in the 
lower panel of Fig. ITT] together with the Rice parameters 
calculated using the surface energy for the { 100 } facet. 
We find that the two sets of Rice parameters (upper and 
lower panels) predict rather similar effects. The situation 
for A is very different: A{100} monotonously increases 
with Ni and decreases with Mn addition. The inconsis¬ 
tency between the Rice parameter and A obtained for 
the { 100 } plane is a consequence of the limited informa¬ 
tion built in the Rice ratio (merely through the surface 
energy) compared to the cleavage stress (involving both 
the Young modulus and the surface energy). We con¬ 
clude that in contrast to the Rice parameter, the average 
A (taking into account both cleavage planes) indicates in¬ 
creased brittleness for Fe-Mn and increased ductility for 
Fe-Ni. For both alloys, the present predictions based on 
A are in line with the observations [l|. 

V. CONCLUSION 

We have investigated the effect of Mn and Ni addi¬ 
tion on the mechanical properties of ferromagnetic bcc 
Fe. Fe-Ni, the lattice parameter increases nearly linearly 
with alloying whereas in Fe-Mn, the lattice parameter in¬ 
creases up to ~ 5 % Mn and then decreases with further 
Mn addition. 

For both Fe-Mn and Fe-Ni alloys, the elastic moduli 
show a nonlinear trends as a function of concentration. 
The surface energy and the unstable stacking fault en- 
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ergy decrease by adding Mn or Ni to Fe. For both planar 
fault energies, Ni shows a stronger effect than Mn. Seg¬ 
regation seems to have a minor effect on the surface and 
USF energies for dilute Fe-Ni and Fe-Mn alloys. Accord¬ 
ing to semi-empirical correlations based on the Poisson 
ratio, the Cauchy pressure, and the Pugh ratio, we have 
found that Mn makes the bcc Fe-based dilute alloy more 
brittle and Ni makes it more ductile. However, the study 
of 7s/7« for the {110} surface and the {110}(111) slip 
systems indicates that both binary alloys become more 
ductile with increasing solute concentration although the 
relative fracture behavior is consistent with the other 
three criteria. This discrepancy between the Rice and 
the Pugh criteria is ascribed to the complex magnetic 
effects around the planar fault, which are however miss¬ 
ing from the bulk parameters entering the Pugh/Cauchy 
criterion. 


Using the theoretical cleavage stress and the shear 
modulus associated with the dominant slip system in 
bcc alloys, we introduce an alternative measure for the 
ductile-brittle behavior. We find that our AjlOO} = 
(Tcz,{100}/G{110}(111) ratio is able to capture the ob¬ 
served alloying effects in the mechanical properties of Fe- 
rich Fe-Ni and Fe-Mn alloys. The superior performance 
of A compared to the Rice parameter lies in the addi¬ 
tional information built in the theoretical cleavage stress 
as compared the surface energy. 
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